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' The energies of 3 H, 3 He, and 4 He ground states, the | and ^ scattering 

q ■ states of 5 He, the ground states of ®He, ®Li, and 6 Be and the 3 + and + ex- 

s : 

cited states of e Li have been accurately calculated with the Green's function 

>' 

: '"j ■ Monte Carlo method using realistic models of two- and three-nucleon inter- 

rS ; 

■ actions. The splitting of the A = 3 isospin P = \ and ^4 = 6 isospin T = 1, 

J n = + multiplets is also studied. The observed energies and radii are gen- 
erally well reproduced, however, some definite differences between theory and 
experiment can be identified. 
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A system of interacting nonrelativistic nucleons is the simplest model of nuclei. Even in 
this simple model exact calculations have been possible only for a limited number of light 
nuclei due to the strong spin-isospin dependence of nuclear forces. For many years, only two- 
nucleon states could be exactly calculated. Next, the Faddeev method was used to study 
the three-nucleon states In the past decade, many advances have become possible due 
to the development of supercomputers. Quantum Monte Carlo methods were used to study 
nuclei with A < 5 HH , the 4 He ground state was calculated with the Faddeev-Yakubovosky 
method f| , and methods using hyperspherical functions were developed to study low-energy 
three- and four-nucleon states |J. In this letter, we report the first realistic six-nucleon (6N) 
quantum Monte Carlo calculations along with updated results for nuclei with A < 5. Until 
now the A = 6 nuclei have been mostly treated as three-body systems with an a and two 
nucleons 0. 

The new Argonne V\s two-nucleon interaction || is used here. It is expressed as a sum 
of four parts: 

V = V 14 + V cib + V csb + V em . (1) 

Its dominant part, t> 14 , contains 14 isoscalar operators as in the old Argonne Vu |J. The 
charge-independence-breaking part, v c n,, has three isotensor terms with operators [3r iz Tj Z — 
Ti-Tj] (g) [1, (Ti-crj, Sy), and includes the effect of the mass difference between charged and 
neutral pions. The v c a, = in isospin T = states, while in T = 1 states v C ib(np) = 
—2v C ib(nn) = —2v C ib(pp). The isovector charge-symmetry-breaking part, v csb , contains the 
operator (r^+r^); it accounts for the difference between nn and pp interactions and vanishes 
in np pairs. The electromagnetic part, v em , contains pp and np Coulomb and magnetic 
interactions in all pairs. The kinetic energy operator associated with this model has isoscalar 
and isovector parts denoted by K and K cs b\ 
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Due to its careful treatment of isospin-symmetry-breaking terms, the new Argonne v±s model 
is well-suited to study the mass differences between the T — |, 3 He — 3 H doublet and the 
T = 1, J 77 = + , 6 He - 6 Li - 6 Be triplet. 

Three- nucleon interactions, V^k, described with the Urbana model [jlO[ are included in the 



nuclear Hamiltonian. These contain a two-pion exchange part with its strength chosen 
to reproduce the observed binding energies of 3 H and 4 He, and a phenomenological spin- 
isospin independent interaction of strength U adjusted to obtain the empirical equilibrium 
density of nuclear matter. In Urbana models IX (VIII) of Vijk, to be used in conjunction 
with the new t> 18 (old these parameters have values A 2n = —0.0293 (—0.028) and 
U = 0.0048 (0.005) MeV. 

The Green's function Monte Carlo (GFMC) calculations ||£|] are carried out with a 
simpler isoscalar Hamiltonian: 

H = K + y £v B (ij)+ E V ijk- (3) 

i<j i<j<k 

The interaction v$(ij) has eight terms, with operators [l,Tj-Tj] <S> [1, (J{ ■ <Jj, Sij, L- S], 
chosen such that it equals the Vu in all S- and P-waves as well as in the 3 Di wave and 
its coupling to the 3 Si. The eigenstates of H are computed, and all other terms, namely 
(t>i4 — V8),v C ib,v cs b,v em , and K cs b, are treated as first order perturbations. The error in the 
binding energy of 3 H due to the perturbative treatment of vu — vg has been estimated by 
Kamada and Glockle to be ~ 0.02 MeV. 



Nuclear states are represented by vector functions ^(R) whose components ip a (R) give 
the amplitudes for the nucleons, at spatial coordinates R = (r 1; r 2 , • • • , r^), to be in spin- 
isospin state \a). The GFMC calculation of the lowest energy state with quantum numbers 
£(= J^yTjTz) starts with an approximate wave function ^(R) determined from a varia- 
tional Monte Carlo (VMC) calculation for that state. It yields values of \I/(t, Rj(r)), where 

^(r,R) = e-^(R), (4) 

at configurations Rj(r) distributed with probability |\L / J ) (R.)\[ / (t, R)|. Here i labels configura- 
tions which number ~ 10, 000. The e~ Hr is considered as a product of many small imaginary 



time steps e~ HAr , with Ar ~ O.OOOlAfeV -1 . Using the small-time limit of the propagator 
e~ H , correct up to order Ar, the f (r + Ar, Rj(r + Ar)) are stochastically estimated from 
the known ^(r, Rj(r)). 

The number of spin-isospin states that contribute to the T = 0(1) 6N ^(R) is 320(576), 
and their propagator e~ HAr is a 320 x 320(576 x 576) complex matrix function. Therefore 
6N GFMC calculations are numerically very intensive. They were performed on the IBM 
SP parallel computer at Argonne National Laboratory with 128 nodes operating at ~ 40 
MFLOPS/node. Propagating 10,000 6N configurations up to r = 0.06MeV _1 requires 
~ 2,000 node hours. In contrast, the A He ground state wave function ^(R) has only 32 
spin-isospin states, and requires ~ 100 node hours for a similar calculation. 

It is well known that there are no A = 5 bound states. The | and | resonances in 
5 He have been studied as n — 4 He scattering and the phase shifts in these partial waves 
have been extracted [|1^]. The GFMC calculations for these states are carried out with an 
external a — n potential that vanishes for \r a — v n \ < 12.5 fm and is infinitely repulsive 
for |r Q — r n | > 12.5 fm. These conditions also imply that the external potential is 0(oo) 
for |r n | < 10 (> 10) fm in the center of mass frame. The calculated \l/ corresponds to a 
scattering solution with a node at |r Q — r n | = 12.5 fm, however, only the interesting part with 
\r a — r n | < 12.5 fm is retained in the calculation. Assuming that the interaction between a 
and n is negligible for |r a — r n | > 12.5 fm the true energy of this state can be obtained from 



the phase shifts [13 



The variational wave functions used in these calculations include spatial and spin-isospin 



two-body and three-body correlations denoted by f c ( r ij),Uij, and £7^ in ref. JT4[ . The 
uncorrelated wave functions for 3 H and A He are as given in [14|, while those for 5 He have 
an additional nucleon in the (p|) or (p|) scattering state. The two extra nucleons in 6N 
states are in optimized superpositions of (p|) 2 , (p|p|), and {p\) 2 states whose radii are 
constrained to reasonable values. Details of the VMC and GFMC calculations will be 
published separately. 
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The transient energy E{r) defined as: 

(gjgffiTj) (tt(j)|ff|tt(§)) , . 

1 j (*„|*(r)> (*(5)|*(I)> ' 1 J 

provides an upper bound which approaches the lowest-energy eigenvalue with quantum num- 
bers £ in the limit r — > oo. The values of E(r), calculated from 10,000 configurations for 
each of the three 6N states and 20, 000 configurations for A He are shown in Fig. [I] along 
with their statistical error estimates. The 4 He E{r) decreases from the variational energy 
E v = E(t = 0) by ~ 1.2 MeV in ~ 0.02MeV , and is almost independent of r thereafter, 
suggesting that the *ff v has < 2% admixture of states with excitation energy > 50 MeV. The 
E{t) of the 6N states, in contrast, decreases by ~ 4 MeV from E v , and does not appear 
to have reached the r — > oo asymptotic value at r = 0.06MeV -1 . The statistical error, 
governed by variance of H^ v (R)/^/ v (Ii) is larger in the 6N E{r). These indicate that the 
6N ty v are not as accurate as the 4 He and that they need to be improved for more 
accurate GFMC calculations. 

The average value of E(t) in the r interval 0.03 to 0.06 MeV~ l is listed as E in Table 
I. The E of 3 H, 4 He, and 5 He can be indentified with their ground state energies, however, 
those of the 6N states can only be regarded as upper bounds. It is difficult to extrapolate 
the 6N E(t) to r — > oo, particularly due to the large statistical errors; nevertheless we have 
attempted it in the following way. Let be the eigenstates with quantum numbers £ and 
energies E^. The ty v contains admixtures of with amplitude fa in addition to the ^ of the 
lowest-energy state. Admixtures with the smallest — E D determine the behavior of E{r) 
at large r. We approximate the contribution of low-lying states with two delta functions at 
E\ and E2 with amplitudes fa and fa. It is experimentally known that the lowest T = 0, 
J n = l + (3 + ) resonances in 6 Li are 5.65(~ 13.6) MeV above the l + (3 + ) bound states fllS" 



Accordingly, we assume that E\ — E Q = 5.65(13.6) MeV for the l + (3 + ) states and also use 
the value 13.6 for T = 1, J 71 " = + . Further assuming that E2 — E a = 30 MeV the calculated 
values of E{r > O.OlMeV" -1 ) are fitted by varying E , fa, and fa. Fortunately, the fits are 
not very sensitive to E2 — E , and are shown in Fig. |l|. The resulting value of E is listed as 



the calculated energy in Table I. The error in the extrapolated E is much larger than that 
in E and less reliable. 

The expectation values of various terms in the Hamiltonian are also listed in table I. 
These are averages over the interval r = 0.03 to O.OGMeV" 1 , calculated using 

<° W> = 2 (*.|*(r» " <6) 

correct up to order |\P(t)) — \^f v )- The dominant (fu) and (i^) have similar values in 4 ife 
and 5 He; while in the 6N states they are closer to the sum of their values in 4 He and 2 H. 
The (Vijk) are similar in 4 ife, 5 He, and 6N states. The (vls) contains both LS and L-S(rj-Tj) 
terms in the NN interaction. They contribute, along with Vijk Jl6|| , to the splitting between 



| and | states of 5 He. However, the calculated splitting of 0.8(3) MeV is much smaller 
than the observed 1.4 MeV. The magnitudes of (vls) and (Vijk) are larger in the 6N + 
and 3 + states than in the 1 + , suggesting that the 1 + has less contribution from the (jof) 2 
configuration than the other two. The underbinding of the 6N + and 3 + states by ~ 1 
MeV is probably related to that of the 5 He | state. In 5 He, only the expectation value 
( v coui) °f the Coulomb interaction between protons is calculated. The other terms in (v em ) 
for 5 He (given in parentheses in table I) are estimated from results for 4 ife and 6 He. The 
( v coui) decreases by ~ 5% from A He to 6 He indicating that the a-cluster expands as we go 
from A = 4 to 6. 

The last three lines of Table I give rms proton and neutron radii. The calculated values 
of R(p) compare well with the those extracted from observed charge radii ||17|| . The exper- 
imental R(p) of 3 He is 1.77 fm, in reasonable agreement with the calculated R(n) of 3 H. 
The Coulomb interaction accounts for most of the isovector 6 Be — 6 He difference (Table II). 
Since this difference is correctly predicted, the R(p) of 6 Be, assumed to be equal to the R{n) 
of 6 if e, appears to have a reasonable value. 

The contributions of v em ,v C ib,v cs b, and K cs b, treated as first order perturbations in this 
work, are responsible for the energy differences within the T = \ and T = 1, J 77 = + 
6N multiplets listed in Table II. The present Hamiltonian explains the isovector energy 
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differences 3 He — 3 H and 6 Be — 6 He fairly well. The three-body calculations show that 
the isovector v cs t, and K cs i } are necessary to obtain the observed 3 He — 3 H difference, in 



agreement with earlier results of Faddeev calculations Unfortunately, the calculated 

value of the isotensor difference |( 6 -Be + 6 He) — 6 Li is much larger than observed. The 
observed difference is essentially explained by the electromagnetic interaction alone. This 
is very puzzling because most of the contribution of the isotensor to this difference 
should be from the relative 1 Sq two-nucleon state in which the difference between pp and 
np phase shifts seems to be well established [19,|20|. The nonperturbative contribution of 
the Coulomb interaction, particularly in e Be, neglected here, may reduce the value of this 
isotensor difference. There could also be some contribution from charge-dependence of the 
two-pion exchange V^. 

In conclusion, we have demonstrated that the GFMC method can be used to accurately 
calculate the energies of the many nuclear states with A < 6 from realistic models of nuclear 
forces. It appears possible to extend these calculations to several A = 7 and 8 states by 
using larger parallel computers such as the IBM SP2. The calculated energies are in good 
agreement with experiment. However, some differences, such as the underestimation of the 
splitting between § and \ states of 5 He are clearly established. We could attempt to 



probe relativistic effects |2i|j2^| , and the spin-isospin dependence of the short-range part of 
the Vijk using these differences. A detailed analysis of the GFMC wave function e~ Hr ^ v is 
in progress to study the structure of the 6N states and improve upon their ty v . 
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FIGURES 

FIG. 1. The top to bottom data sets show the E(t) for 6 He, 6 Li(3 + ), 4 He, and 6 Li(l+) states, 
along with the fits described in the text. Constants of 2.5, 1.5, and 0.5 MeV have been added to 
the E(t) of 6 He, 6 Li(3 + ), and A He, respectively, for clarity. 
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TABLES 



TABLE I. Calculated energies and radii in MeV and fm. 



Nucleus (J) 


2 H(l) 




4 He(0) 


5 tfe(§) 




6 He(0) 


6 Li(l) 


6 Li(3) 


E(Expt.) 


-2.22 


-8.48 


-28.3 


-27.2 


-25.8 


-29.3 


-32.0 


-29.8 


E(Calc.) 


-2.22 


-8.47(2) 


-28.3(1) 


-26.5(2) 


-25.7(2) 


-28.2(8) 


-32.4(9) 


-28.9(6) 


E 


-2.22 


-8.47(2) 


-28.3(1) 


-26.5(2) 


-25.7(2) 


-27.3(4) 


-31.1(4) 


-28.2(3) 


(K) 


19.9 


50.(1) 


118.(1) 


122.(2) 


117.(2) 


146.(4) 


143.(3) 


138.(3) 


(vu) 


-22.1 


-58.(1) 


-142.(1) 


-145.(2) 


-140.(2) 


-172.(4) 


-173.(3) 


-165.(3) 


(V ijk ) 


0. 


-1.20(3) 


-6.5(3) 


-7.0(4) 


-6.4(4) 


-7.0(7) 


-6.2(6) 


-6.9(5) 


(VLS) 


-0.08 


-0.20(5) 


-0.4(1) 


-1.2(1) 


-0.4(1) 


-2.7(3) 


-1.5(5) 


-3.0(4) 


(Vem) 


0.018 


0.039(1) 


0.879(5) 


(0.86) 


(0.86) 


0.87(1) 


1.71(2) 


1.72(2) 


&) 


0. 


0. 


0.761(2) 


0.745(3) 


0.751(3) 


0.724(8) 


1.55(2) 


1.57(2) 


R(n) 


1.967 


1.72 


1.42(1) 


3.02(3) 


3.57(3) 


2.62(1) 


2.41(5) 


2.46(7) 


R(p) 


1.967 


1.58 


1.42(1) 


1.84(2) 


1.99(2) 


1.89(6) 


2.41(5) 


2.46(7) 


R(p)(Expt.) 


1.953 


1.61 


1.47 








2.43 





TABLE II. Energy differences in MeV within isospin multiplets 

3 He- 3 H 6 Be- 6 He ±( 6 5e + 6 He) - 6 Li 



(v em ) 0.669 2.33 0.33 

(K csb ) 0.014 0.036 

(v csb ) 0.066 0.116 

(va b ) 0.28 

A(Calc.) 0.749(1) 2.5(1) 0.6(1) 

A(Expt.) 0.764 2.35 0.34 
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• • 6Li 1+ 




0.00 0.01 0.02 0.03 0.04 0.05 0.06 



x (1/MeV) 



